Generic phase diagram of active polar films 
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We study theoretically the phase diagram of compressible active polar gels such as the actin 
network of eukaryotic cells. Using generalized hydrodynamics equations, we perform a linear stability 
analysis of the uniform states in the case of an infinite bidimensional active gel to obtain the 
dynamic phase diagram of active polar films. We predict in particular modulated flowing phases, 
and a macroscopic phase separation at high activity. This qualitatively accounts for experimental 
observations of various active systems, such as acto-myosin gels, microtubules and kinesins in vitro 
solutions, or swimming bacterial colonies. 
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Active materials are a challenging class of systems 
driven out of equilibrium by an internal or an exter- 
nal energy source. Examples of active systems are self- 
propelled particle assemblies in bacterial colonies [l], , 
or the membrane and the cytoskeleton of eukaryotic 
cells Q . The cell cytoskeleton is a network of long fila- 
ments made by protein assembly, interacting with other 
proteins Q which can, among other things, crosslink or 
cap the filaments. Motor proteins, myosins, kinesins or 
dyneins use the chemical energy of Adenosinctriphos- 
phate (ATP) hydrolysis to "walk" along the filaments, 
and exert stresses that deform the network [5j, leading 
to an active behavior. The active character of the cy- 
toskeleton plays a major role in most cell functions such 
as intracellular transport, motility and cell division. 

The cell cytoskeleton has a rich and complex dynami- 
cal behavior @, d, 0, B @] ■ Self-organized patterns, such 
as asters, vortices, and rotating spirals, microscopic and 
macroscopic phase separations ("superprecipitation" 
) have been observed as a function of motor and ATP 
concentrations in a thin film [Bj|. This two-dimensional 
geometry gives for example a good description of the thin 
lamellipodium of a cell spreading or moving on a sub- 
strate. Some of these structures have been recently ex- 
plained theoretically [ll|, ll2j , but a full phase diagram of 
active polar films is still missing. 

The passive visco-elastic properties of the cytoskele- 
ton are similar to that of a physical gel made of the 
cross-linked semi-flexible filaments. Recently, Kruse et 
al. 12| have proposed a generalized hydrodynamic the- 



ory to describe macroscopically the active character of 
incompressible polar gels, based on conservation laws 
and symmetry considerations. In this letter, we use the 
generic model of Ref.[12|] to study the stability of com- 
pressible active polar films. We perform a linear stability 
analysis of the uniform states in the case of an infinite 
two-dimensional geometry, and obtain the dynamic phase 
diagram. Our results qualitatively account for the exper- 
imental observations on various active systems, such as 



acto-myosin gels, microtubules and kinesin solutions in 
vitro or swimming bacterial colonies. We choose the ex- 
ample of infinite acto-myosin films that we consider as 
two-dimensional, and non interacting with the environ- 
ment. Such an example could be realized by a freely 
suspended film such as those studied for liquid crystals. 
We consider the long wavelength limit, and all variables 
are implicitly averaged over the film thickness. The actin 
network has a local macroscopic polarization given by a 
unitary vector p = (cos 8, sin 6). This describes degen- 
erate parallel boundary conditions at the film interface, 
but also equivalently normal boundary conditions with a 
splayed state. Since the film thickness can vary, the two- 
dimensional actin network is compressible, even though 
the three-dimensional material is not. The gel has a 
weakly fluctuating density c(r) = cq + p(r). The aver- 
age density cq can be set to 1 by rescaling the various 
coefficients of the free energy given in EqQ] below. As a 
first approximation we neglect the interactions with the 
solvent and use a one fluid model. 

The free energy of the gel, up to quadratic order, cou- 
ples the polarization p to the density fluctuation p : 
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The splay and bend elastic moduli are assumed to be 
equal (K\ = Ks = K) for the sake of simplicity. The La- 
grange multiplier h\\ enforces the constraint p 2 — 1. The 
spontaneous splay term leads to boundary terms that are 
irrelevant in the infinite system size limit (l3j . The vari- 
ation of the free energy with density is characterized by 
the positive compressibility a and by the positive coeffi- 
cient (3 associated to the density fluctuation correlation 
length; w is a coupling constant between density fluctu- 
ations and splay. The molecular field, conjugate to the 
polarization is h a = —5F/5p a with coordinates (hu,h±) 
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parallel and perpendicular to the polarization. The free 
energy is similar to that of a ferro-electric nematic liquid 
crystal when the order parameter n does not have a fixed 



length [lj,115l) 



The gel motion is described by the two-dimensional 
velocity field v. The strain rate tensor is u a p — (d a vp + 
d/3V a )/2, and the vorticity tensor uj a p — (d a vp — dpv a )/2. 
The conservation equation of the gel is written as: dtp + 
d a (l + p)v a = 0. 

The gel is driven out of equilibrium by continuous and 
homogeneous input of energy, characterized by the chem- 
ical potential difference Ap between ATP and its hydrol- 
ysis products, which we assume to be constant. 

The gel dynamics is described by the linear hydrody- 
namic equations for active polar gels of Ref.[l2|. The 
constitutive equations for the mechanical deviatory stress 
tensor a a p and the rate of change of the polarization 
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when the gel behaves as a viscous liquid: 
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We neglect here the geometric non-linearities introduced 
in (l2j . The rotational viscosity 7 and the coupling con- 
stants between flow and polarization v, D are standard 
liquid crystal parameters [161 ]. The active contributions 
to the mechanical stress and to the rate of change of the 
polarization are proportional to Ap and are character- 
ized by the coefficients £, £ and A. This set of constitu- 
tive equations is completed at low Reynolds numbers by 
the force balance: d a (a a p — HS a p) — 0. Locally, there are 
two forces acting on the gel, the deviatory stress tensor 
a a f3 and the pressure II = ^ = toV • p + ap — (3Ap. 

It is observed experimentally [171 ] , and predicted by 1- 
dimensional models |6j that the overall effect of myosin II 
motors on actin solutions is contractile. This corresponds 
to negative values of both £ and The molecular motors 
also have an effect on the rate of change of the polariza- 
tion described by A; if A > the polarization is enhanced. 
This seems to be observed experimentally (zipping effect 
in Ref.[l8|). We consider here only ordered polar phases 
with a unitary polarization vector p, and we ignore this 
active coupling (setting A = 0) for simplicity (l9j . 

In order to discuss the accessible steady states of an in- 
finite active gel film, we first study the states of uniform 
polarization and velocity gradient. These asymptotic 
states are non-equilibrium states, and cannot therefore 
be obtained by minimizing a free energy functional. In a 
two-dimensional geometry, the hydrodynamic equations 



of motion give 11 scalar equations for the 11 indepen- 
dent variables x = 6, p, h a , v a , <7 a p, it. These equations 
of motion have two types of homogeneous steady states: 
a static state where the velocity gradient u a p vanishes 
and the polarization is uniform and oriented in a direc- 
tion #0 ; and a flowing state with a finite velocity gradient 
and a uniform polarization. In two dimensions, the most 
general steady flow with a constant velocity gradient is 
a superposition of two simple shear flows in two perpen- 
dicular directions x and y. The only non vanishing com- 
ponents of the velocity gradient are d x v y and d y v x . We 
consider here for simplicity only one component shear 
flows for which d y v x — 0. The velocity is then along 
the y direction and the polarization angle 9 V is such that 
cos(26*„) = \jv (we assume \v\ > 1). This flowing state 
is the analog, for an infinite compressible gel, of the flow- 
ing state obtained in a confined geometry in 13| and to 
the rotating spirals in a cylindrical geometry described in 
12]. It confirms the possibility of obtaining spontaneous 
flows in polar active materials. 

The stability of the two homogeneous states is studied 
by introducing a small perturbation at point r, at time 
t with a wavevector k and a growth rate s: x = xo + 
Xi exp[st + zkr], where Xo is a steady state solution. The 
equation for the perturbation xi can formally be written 
in matrix form Mxi = and the possible growth rates 
of the perturbation are determined from the equation 
det(M) = 0. This is a quadratic equation in s with 
two roots denoted by s + , s~ with Re(s + ) > Re(s~). The 
sign of Re(s + ) gives the stability limit of the steady state 
homogeneous phase. 

We are only able to give a complete discussion of the 
stability of the uniform steady states with respect to a pe- 
riodic perturbation in the quasi 1-dimensional case where 
we do not allow for a y dependence (k y = 0). The treat- 
ment of the fully general 2-dimcnsional problem requires 
numerical work. However most of the physics can be 
extracted from the quasi- 1-dimensional case, which we 
present hereafter; the numerical study allows us to ex- 
tend our conclusions to the general 2-dimensional case. 

We first analyze the stability of the static state, fix- 
ing all the parameters of the active gel introduced above, 
except for the Franck constant K and the active stress 
CAp. As there is no flow in the system, there is one 
single direction, given by the polarization, which we 
choose as the y axis (#0 — 7r /2)- The numerical study 
of the 2-dimensional problem reveals that the maximum 
of Re[s + (k x , ky)] lies on the axis k y = 0, for any value 
CAp < (Ap*(K), with C,Ap*(K) > 0. In this regime, 
it is therefore sufficient to consider one-dimensional per- 
turbations in the x direction. The growth rates are then 
solutions of: 
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2rj-fas 2 + s [(2r}/3 + Kab-/)k 2 + 2rja + C,A^a{v + 1)] + bK/3k 4 + [b(Ka - w 2 ) + (A(j,/3(v + 1)] k 2 + (A^a(u + 1) = 



where we have redefined k x = k and where a and 
b are dimensionless functions of the parameters: a = 
2^ + v(v + v) and b = 2^ + fi±^. It is useful to ana- 
lyze first the passive case where (Apt = 0. If K > w 2 /a, 
Re(s + ) is maximum and negative for k = and the active 
gel is stable. If K < w 2 /a, one can check that Re(s + ) 
is maximum and positive when k = k c ^ 0. The uni- 
formly ordered phase is locally unstable with respect to 
a finite wavelength longitudinal mode transverse to the 
ordering direction. It follows that in this region of the 
phase diagram there necessarily exists an ordered modu- 
lated phase characterized by a periodic polarization. The 
precise symmetry of this passive phase, (striped, hexag- 
onal or cubic phase), has been discussed in [3, EH for 
the case of ferroelectric nematic liquid crystals. We now 
consider the effect of the activity on this instability. For 
(A/j, < 0, the static state is unstable at zero wave vec- 
tor (both growth rates are real and s + s~ < 0], with 
respect to the flowing state, as found in Ref.[13j. For 
CAfi > 0, the uniform state is unstable with respect 



to a finite wavelength longitudinal mode, transverse to 
the polarization direction when the Franck constant is 
small enough K < K c (C,A^l). The function K c ((A/j,) 
can be analytically calculated from EqH] in this regime 
and is plotted in Fig[TJ For CA^ > (Afx*(K) > 0, a 
2-dimensional numerical analysis is required and reveals 
that for K < K c (CAfi) (evaluated by numerically com- 
puting the sign of Re(s + )), there are 2 independent most 
unstable wavevectors k c i,and k C 2. 

We now discuss the flowing state ((Afi < 0), which 
we believe to be relevant for cytoskeleton dynamics. (If 
CAfi > 0, the flowing state is unstable at any wave 
vector). The numerical study of a 2-dimensional per- 
turbation shows that the maximum of the growth rate 
Re[s + (fc x , k y )] is obtained for a wavevector k perpendic- 
ular to the polarization, for any value QAfi < 0. It is 
therefore sufficient to consider an effective 1-dimensional 
problem where spatial variations are allowed only along 
the unstable direction orthogonal to p. The growth rates 
are in this case solutions of : 



2rjjas 2 + s [(2r](3 + Kab^k 2 + 2r)a'] + bK/3k 4 + b(Ka - w 2 )k 2 + ikdw(A[j, - ^Ajia^ + 1) = (5) 



where a' is an effective compressibility, and d a dimen- 
sionless function of the parameters: d = 2^ + (^+l) /2 + D 
and a' = a + d-yCA^v + l)/2. 

For K < K c ((Afi) (evaluated by numerically comput- 
ing the sign of Re(s + ) from Eq[5]), the uniform state is 
unstable with respect to a finite wavelength mode. The 
value of the critical Franck constant K c diverges for an 
active stress (A/jL c such that a'(CA/i c ) = 0. At lower 
values of the active stress, the flowing state is always 
unstable. The most unstable wave vector |k c | decreases 
with £A^t and vanishes if the active stress is smaller than 
a critical value C,A^, p (K). Note that in this regime of 
negative active stress, Im(s + ) ^ and the instability is 
oscillatory, as opposed to the regime £A/i > 0. 

The results of the stability analysis are summarized in 
the phase diagram of FigfT] A full analysis of the non- 
linearities is necessary to predict the symmetry of the 
dynamic equilibrium states. This would require an ex- 
haustive study of all the possible quadratic and cubic 
terms in the perturbations in the equations of motion, 
which seems out of reach. We can nevertheless infer qual- 
itatively the structure of the different phases. 



In regions B _ and B + of the phase diagram, the in- 
stability occurs at finite wavevector. It is similar to the 
instability of the passive system described in 1J, |l5( : it 
favors splay in the system, which can appear either in 
a striped or "lattice" (cubic, hexagonal, or triangular) 
phase. Because of the outbreak of polarization gradi- 
ents, the strain rate tensor cannot vanish and the sys- 
tem flows. The continuity of the velocity field imposes 
extra constraints and does not allow for the formation 
of domain walls. This, together with the results of [l2| 
suggests that in the lattice phases, each elementary cell 
contains a spiral-like structure with a rotating flow. This 
structure is not compatible with the hexagonal symme- 
try : the only possible lattice phases have therefore cubic 
or triangular symmetries sketched on Fig[TJ 

In region C~, the system is unstable for k = and 
thus its behavior depends on boundary conditions. We 
expect a macroscopic phase separation since the effective 
compressibility a' becomes negative. In region C + , the 
existence of 2 independent unstable modes suggests that 
an oblique phase prevails, even if non linear terms could 
a priori select only one of the unstable modes. 
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FIG. 1: Dynamic phase diagram of an active polar film. 
The transition line between regions B + and C + is the line 
£An*(K); the transition between regions B~ , B + and C + 
and regions A~ and A + is the line K c (^Ajj,); the transition 
between regions B~ and C~ is the line A/j, p (K) 



therefore polar; in a coarse-grained picture, a colony can 
be viewed as a viscoelastic gel; last, they consume chem- 
ical energy (O2) and hence are active. The alternative 
model of [l| also predicts an intrinsic flow instability for 
self propelled particle assemblies. This instability differs 
from the one discussed here in that it has no thresh- 
old. The main differences with our description are first 
that our model in this one fluid version cannot impose a 
spontaneous velocity of the gel with respect to the back- 
ground fluid; second, our model describes a compressible 
gel, this feature being crucial to trigger the instabilities. 
A thorough comparison between the two theories needs 
a two fluids description and is underway. 

We are grateful to F. Amblard (Institut Curie), D. R. 
Nelson (Harvard University), D. Riveline (Grenoble) and 
S. Ramaswamy (Bangalore) for very useful discussions. 



Our analysis has been presented in terms of the actin 
cytoskeleton, but it is general enough to be applied to any 
visco-elastic, polar, and active material. We now briefly 
discuss two other examples. One should keep in mind 
that our model does not include noise, be it thermal or 
intrinsic (due to the stochastic activity of motors, actin 
polymerization...) and that the predicted ordered phases 
could be disordered in the presence of noise. 

Systematic quantitative data describing the phases 
of in vitro acto-myosin solutions as a function of the 
contractility (associated to changes either in myosin II 
and/or ATP concentrations) and therefore of £A/z do not 
seem to be available yet. However, both disordered vor- 
tex flowing phases [111, G3 , which could correspond to re- 
gion B~, and phase separations (superprecipitation, [loj ] . 
corresponding to region C~) are observed in experiments 
or obtained by numerical simulations. Microtubules and 
kinesins solutions also belong to the class of materials 
described by our model. They are viscous, polar, and ac- 
tive. The phase diagram obtained in this paper qualita- 
tively agrees with the phases observed in vitro for micro- 
tubules and kinesin motors in two dimensions when 
the motor concentration is increased, a flowing phase of 
spirals appears (B); for higher motor concentrations, the 
spirals become growing asters which progressively sepa- 
rate (C~); for even higher concentrations microtubules 
bundles are formed (C~). 

A similar behavior is observed in bacterial colonies : 
experiments on 2-dimensional colonies of B. SubtilisQ 
show a bacterial flow composed of rotating whirls of 
swimming bacteria. This bacterial turbulence could be 
described by a disordered version of the cubic phase (B+ , 
C + or B~) predicted by our model. Indeed these bacte- 
ria are rod-like shaped, oriented by their flagella, and 
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